Exact analysis on the singularity of the joint density of states and its relationship to 

the quasiparticle interference 
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Singularities of the joint density of states (JDOS) and Fourier-transformed local density of states 
(FT-LDOS) correspond to the hot spots in quasiparticle interference patterns. In this paper the sin- 
gularity of JDOS is analyzed exactly, with three types of singularities being classified. In particular, 
the third type of singularities are found exactly to be envelopes of the contours of constant energy. 
Remarkably, we show that JDOS and FT-LDOS have the same singular points. Approaching to 
the singular points, both quantities diverge complementarily in an inverse-square-root manner if the 
joint curvature is nonzero. The relative magnitude of divergence is governed by the joint curvature 
as well as the product of the quasiparticle velocities. If the joint curvature of certain singularity is 
zero, the divergence has a higher order than —1/2. 
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INTRODUCTION 

Fifteen yeas ago Fourier-transformed (FT) scanning 
tunneling microscopy/spectroscopy (STM/STS) first 
demonstrated its potential as a supplementary technique 
to probing the surface band structure of metals jlj . It is 
astonishing that the intrinsic momentum-space electronic 
structure can be reflected by imperfections and obtained 
by a real-space measurement. Since then the FT-STS 
technique has been widely applied to study the Fermi 
surface and energy gap structure in <i-wave cuprates0- 
4] and iron-based superconductors [jj H as well as the 
surface bands of topological insulators [7H9| by imaging 
the modulation of local density of states (LDOS). The 
modulation is caused by the interference of quasiparti- 
cles scattered elastically off imperfection on the material 
surface. If the quantum interference is between quasi- 
particles scattered from to kj states on the contours 
of constant energy (CCE), then the LDOS modulation 
with wave vectors q = kf — k^ appears. Based on the 
quantum transition theory it was argued |2j,|3j that the in- 
tensity of modulation with wave vector q is proportional 
to the joint density of states (JDOS) of the initial and 
final states, i.e. |V25k i | -1 |V25k / | ■ In <i-wave cuprates 
due to the strong anisotropy of the energy dispersion of 
Bogoliubov quasiparticles, the CCE are banana shaped 
with large density of states on their tips. Therefore the 
scattering of quasiparticles from tip to tip on the CCE 
are strengthened according to the JDOS picture, which 
interprets the main findings about the quasiparticle in- 
terference (QPI) pattern in high-T c superconductors. 

The relevance of the JDOS picture to FT-STS has 
been suggest ed [l(| theoretically and checked experimen- 
tally [TkHL Numerical H H 01 

and approximately 



pattern actually corresponds to the FT-LDOS. JDOS 
and FT-LDOS are obviously different physical quantities 
and the matching between them is rather curious. 



analytical 3, 15 1 works show that JDOS surely can re- 
peat the QPI pattern observed by FT-STS and match 
the observed octet vectors quite well. However, the main 
inconsistency of the JDOS interpretation is that the QPI 



Under the Dirac cone approximation, Refs. [14| and 
showed analytically that JDOS and FT-LDOS di- 
verge on the same singular boundary on which the octet 
vector is located. These findings gave a clue on the va- 
lidity of JDOS picture in explaining the QPI pattern. 
That is, the hot spots in the QPI pattern observed by 
FT-STS possibly correspond to the singularities where 
both JDOS and FT-LDOS strongly diverge. However 
these works are limited in that (1) it is only applicable 
to (i-wave superconductors in low energy limit where the 
Dirac cone approximation is valid, (2) of all the octet 
vectors only the analytic result for qi can be given and 
(3) one need to assume that the dominant impurities are 
nonmagnetic to highlight qi from other points on the 
boundary. If the scatterer is of magnetic or pair-breaking 
type, qi can even become a dark spot according to the 
results under the Dirac cone approximation (IB) . Actu- 
ally there are various types of scatterers on the surface 
of high-T c cuprates and even a single impurity can scat- 
ter Bogoliubov quasiparticles in a complicated and com- 
posite manner. (4) The JDOS arguments seem to only 
applicable to superconductors with strong anisotropy in 
the energy dispersion. 

The motivation of this paper is (1) to give a concrete 
theoretical foundation for the JDOS picture by explor- 
ing the singularities of JDOS and FT-LDOS exactly, and 
(2) to make the JDOS picture also applicable to more 
general systems. The main results of this paper are as 
follows. (1) The singularities connect two momenta at 
which the quasiparticle group velocities are (anti)parallel 
to each other. (2) For any given quasiparticle energy, the 
singularities form closed curves which are envelopes of 
the one-parameter families of CCE curves. (3) JDOS ex- 
hibits power- law divergent behavior with exponent —1/2 
if the singular point is approached from one side. The co- 
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efficient of the divergent term is governed by the inverse 
of joint curvature as well as the quasiparticle velocities. 
(4) If the joint curvature is zero the power-law exponent 
is higher. (5) The FT-LDOS shares the same singulari- 
ties of JDOS and the divergent behavior is also the same 
except that it diverge if the singular point is approached 
from the other side. 

The main results of our paper are given by Eqs. ([2]). 
(O, (|T2"j) and (fT5|) . Some special cases are discussed to 
show how to make use of these results. In the following 
section we first study the JDOS which is easier than FT- 
LDOS to reveal their singularities. 



SINGULARITIES OF JDOS 

The JDOS, for studying the QPI pattern on the surface 
of metals, superconductors, topological insulators and so 
on, is defined as, 

J (p, uj) = J F kik+V 5(uj - E k+p )S(uj - E^)dk x dk y , (1) 

where E^ = E{k x ,k y ) is the energy dispersion of quasi- 
particles on the surface. -Fk,k+p is a regular function and 
takes into account the effects of coherent factors 15l-ll7 
if the quasiparticle is a superposition of an electron and 
a hole, or phase factors if the quasiparticle is a superpo- 
sition of a spin-up and spin-down electrons. 

We first change the variables from (k x ,k y ) to (E,E') 
with E = Ek, E' = -Ek+p and the double integral in 
Eq. ([T]) is transformed step by step to, 



f 5(uu - E')8(uj - E) , , , 
J(p,w)= / Fk,k+p V — „ ; L dEdE' 



k n + P 



d(k x k y ) 



3{k x k"y 



(2) 



k=k 



E 



-Pkp.ko+p 
Vk X V ko+p | 



where Vk = Vk£k the group velocity of the quasiparticle 
and the sum is over all ko's in the real solution set §(p, uj) 
of the following two equations, 



£k + P = uj. 



(3) 
(4) 



It can be seen that Eq. ([3]) and ^ determine two sets of 
CCE with the latter is just the former translated by —p. 
Geometrically ko 's are all the intersection points of the 
two CCE. As shown in Eq. the point (p, uj) in the 
3D momentum-energy space is a singularity of JDOS if 
there exists at least one k € §(p,w) satisfies [18J, 



Vk x v k +p 



0. 



(5) 



According to Eq. ©, ® and ©, we can clas- 
sify three types of singularities of JDOS. The first one 
is the von-Hove singularities (vHS) which are also the 
well-known singularities of the density of states (DOS), 
i.e. D(uj) = J 8(uj - E k )dk = J Ek=uj |vk| -1 dZ. The vHS 
of D{uj) occur at uj ci for which the CCE crossing the 
critical points k cr with Vk cr = 0. Similarly as indicated 
by Eq. © the vHS of J(p,u>) also occur at the same 
uj ci of the DOS if k or k + p coincide with the critical 
points. Corresponding to uj ci the singularities of JDOS 
in the momentum space is given by the implicit function 
^ko-ip = uj ct . 

The second type is called perfect-nesting singularities, 
which means that there are uncountably infinite num- 
ber of ko satisfying Eq. ([3]) and (j4|) for certain vectors p 
and uj. Geometrically this case means that the two CCE 
have finite segment in common, and thereby Eq. (J5]) is 
naturally satisfied on the common arc since Vk is normal 
to the CCE. The perfect-nesting type of singularities is 
quite special and occurs accidentally in real materials ex- 
cept for the trivial self- nesting case with p = 0. Actually 
this type is the limiting case of the third type, i.e. the 
nearly-nesting type, of singularities. In the following, we 
will elaborate on the third type of singularities, which 
exists inevitably in real systems and therefore is more 
general and physically more relevant. 

If (p,w) is a third type of singularities of JDOS, then 
the coupled Eqs. ([3]), (j4]) and (O can be satisfied un- 
der the conditions: (i) the two CCE do not cross any 
critical point, (ii) the number of ko, i.e. the number of 
crossing points of the two CCE, is finite. Geometrically 
this means that the two CCE are tangent to, or in other 
words touching with, each other at certain point ko as 
indicated by Eq. ([5]). Furthermore the tangent point is 
actually the degenerate root of Eqs. (j3|) and (j4}, whose 
multiplicity is an integer greater than 2. This multiplicity 
can be assigned to every singular point (jp,uj) to charac- 
terize the degree of nesting between the two CCE at ko 
and ko + p, and we denote it by M(p,w). The above- 
mentioned second type, i.e. the perfect-nesting type, of 
singularities corresponds to M(p,uj) = oo. The JDOS 
shows power-law divergence when any singularity (p,u>) 
is approached and the exponent is M(p,w) _1 — 1 as we 
reveal in the next section for the M = 2 cases. 

After canceling ko of the coupled equations ([3]), (UJ 
and ([5]) formally, one can in principle obtain equation for 
these non-trivial singularities of J(p, uj). For any given 
uj these singularities can be expressed as implicit curves 
T™(p) = with to the curve label, each of which is the 
envelope of a single-parameter family of curves (See Ap- 
pendix A for detail). Actually the CCE determined by 
Eq. (J3j> are generally composed of ./V disconnected sim- 
ple closed curves labeled by n = 1, 2, ..N (e.g. in d-wave 
cuprates the CCE consist of N = 4 banana-shaped curves 
in certain energy range). The n-th curve of the CCE can 
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be represented by a parameter equation. For every pa- 
rameter of the n-th curve, Eq. (fJJ determines N families 
of closed curves in the 2D p plane. Therefore there are 
totally N x N families of curves. At last Eq. ([5]) de- 
termines the envelope function of every family of curves, 
i.e. T™(p) = 0. Since the CCE are closed curves, the sin- 
gularities T™(p) = are also composed of closed curves. 

These closed envelope curves then separate the p plane 
into different regions and serve as their boundaries. For 
any p in these regions the intersection points of Eq. §3§ 
and Eq. ((4]) (i.e. the real roots of these two coupled equa- 
tions) are distinct, and the number of intersection points 
are even (0, 2, 4, ... ). Crossing the boundary (singulari- 
ties) from one region to another, the number of distinct 
intersection points changes by 2, because on the bound- 
aries there are degenerate roots (touching point) . 

Eqs. ([2]) and ([5]) also show the significant choosing ef- 
fect of the singularities. Namely on the singular point p 
only the factors i"ko,k +p is highlighted. Therefore the 
variation of -Fk,k+p can be reflected in the JDOS data. 
This effect has been utilized to detect the gap symme- 
try of d-wave cuprates and iron-based superconductors 
via FT-STS 0, As we show in the following sec- 

tion, FT-LDOS and JDOS share the same singularities. 
Therefore FT-LDOS can bear the same choosing effect 
as JDOS. 



DIVERGENT BEHAVIOR OF JDOS AND 
FT-LDOS 

It is easy to propose the idea of singularities from the 
aspect of JDOS as shown above. However the singular 
behavior at the vicinity of the singularities is difficult to 
obtain from Eqs. ©, ©, (gj and ©. This problem 
and the relation between JDOS and FT-LDOS will be 
tackled in this section from another view. We suggest 
that the hot spots in the QPI pattern as observed by FT- 
STS are singularities of FT-LDOS, which can be captured 
by " 
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S(p,u) 



-Fk,k+p(^) 

U — £V + r> 



S(oj — Ek)dk x dk y , (6) 



As in Eq. (fTJ) -Fk,k+p(w) is responsible for the effects from 
nature of two or more components of quasiparticles and 
does not give rise to the singular behavior of FT-LDOS. 
In fact these singularities of FT-LDOS and JDOS are 
embedded in their elementary mathematical structures, 
i.e. the rational and S functions from real and imaginary 
parts of single-particle Green's function. Therefore we 
will temporarily omit -Fk,k+p(w) and concentrate on the 
singularities intrinsically originating from the energy dis- 
persion of the quasiparticles. i<k,k+p(w) will appear in 
the final results of FT-LDOS as a numerator, just like it 
does in Eq. ©. 



We first combine FT-LDOS with JDOS to form a com- 
plex response function as we did before [13, EH , 



R(p,u>) = S(p,u>) + inJ(p,oj) 
1 



— S(uj - E k )dk x dk y , 



(7) 



with the help of which divergence behavior of and the 
relation between JDOS and FT-LDOS can be simultane- 
ously obtained. The 2D integral is transformed to line 
integral along the CCE, 



R(p,u) 



1 



1 



E k =uj 



E k+P + iO+ |V k £ k 



■dL 



(8) 



For any given ui the CCE are generally composed of N 
simple closed curves, each of which can be labeled by an 
integer n = 1,2, ...,N. As we assumed in the previous 
section the CCE do not cross any critical point. There- 
fore each closed curve can be parameterized by a peri- 
odic function. Here we use the arc-length parametriza- 
tion followed by re-scaling of the period to 2tt. There- 
fore the explicit function of the curve can be written as 
k™M = (/"(*)> 9u{t))> where n is the label of the closed 
curve, /, g are regular and periodic functions with period 
2tt. dl n (t) — -Tf-dt' with P" the perimeter of the n-th 
closed curve. After the parametrization R(jp, uS) can be 
expressed as sum of line integral along every single closed 
curve, 



N pn f 



f 



dt 



Ex 



k£(*)+P 



(9) 



where v™(t) = | Vk-E-k|k=k»(t)' In the following derivation 
we will omit the index n for brevity. 

By observing Eq. (JSJ one can find that the integral 
might be singular if the denominator of the integrand 
becomes zero at certain to, that is 



Ex 



M*o)+p 







(10) 



However if to is a distinct root, i.e its multiplicity is 1, 
it contributes a finite value to the integral. Actually 
Eq. pH)) is the parameterized version of Eq. ^ coupled 
with Eq. (|3|) . Therefore the distinct root is equivalent to 
distinct intersection points of Eq. ([3]) and (HJ) . The inte- 
gral is singular only if one of the real roots of Eq. dTU|) is 
degenerate, i.e. 



Ex 



dEx 



Mt )+p 



k„(t)+p 



dt 



= 0. 



(11) 



In fact, Eq. (fTTj) can be easily proved to be equivalent to 
Eq. ©. Therefore (p,w) is a singularity if Eq. ()10p and 
(fTTj) can be satisfied by certain real to. 

For any singularity, R(p, uS) is infinite since the inte- 
grand has a second-order pole to (see Appendix B for 
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detail). What we are interested in is how it diverges to the behavior of R(p + S,ui) with \S\ ~ 0. The result is, 

infinity when approaching to the singularities, namely (see Appendix C for detail) 

I 

D , , x s 1 1 j sgn(n ko+p ■ eko+p,k ), <5 • £k 0+P ,k > 

|vk ||v ko +p| a • eko+p,kol 7 h, 5 • £>k +p,ko < o 



where we define the joint curvature, 

£k',k = Kk'nk' - K k n k . (13) 

Here Kk is the curvature of the CCE at k, whose defini- 
tion is given in Appendix D. rik = Vk/|vk| is the direc- 
tion vector of the group velocity, ko is the real solution 
of Eqs. ©, and ©. 

From Eq. (fT2")) we have the following results: (1) The 
so-called JDOS argument for d-wave superconductors 
is confirmed that the QPI pattern is proportional to 
(l v i|l v /l) _1 although it did not point out that the two 
velocities should be (anti)parallel as we show here. (2) 
There is a more important factor governing the intensity 
of divergence, i.e. the joint curvature fPki.k/. (3) As real 
and imaginary parts of R(p, w), FT-LDOS and JDOS 
share the same singularities and diverge complementar- 
ily in a inverse-square-root manner if any singularity is 
approached, unless the joint curvature is zero. (4) After 
a similar but tedious derivation as we did in Appendix B 
and C, we find that if the joint curvature of certain singu- 
larity is zero, the exponent of the power-law divergence is 
greater than —2/3. In fact the exponent is related to the 
multiplicity M of the touching point by M^ 1 — 1, with 
M = 2, 3, oo. Larger value of M means higher order 
of divergence and better nesting of the two CCE at ko 
and ko + p. And the infinite limit case corresponds to 
the perfect-nesting type of singularities. 

DISCUSSION FOR SOME CASES 

Near its singularities, the FT-LDOS 5(p,w) has 
stronger intensities and are reasonably easier to be ob- 
served by the FT-STS technique than regular ones. Fur- 
thermore higher-order singularities, if they exist, are 
more robust than lower-order ones in the presence of 
various destruction effects on QPI, such as random im- 
purity distribution, multiple impurity scattering, as well 
as short coherent length/time of quasiparticles. In the 
following we will show how to make use of Eq. (TT2"j) to 
find the higher-order singularities. 

If the CCE consist of one single curve whose curvature 
does not change sign, then the singularities of S(p,w) 
are all —1/2 type although the magnitudes may be dif- 
ferent depending on the detailed energy dispersion. How- 
ever if the curvature of the CCE change sign, there must 




FIG. 1. Certain CCE for a toy model. By observation, we 
can see that the curvature of this CCE changes sign, pi, 2, 3 
are singularities of S(p, cj) because they connect points with 
their group velocities (anti)parallel. According to Eq. (|13|) . 
only the joint curvature of p3 is zero which leads to the fact 
that the order of pa is higher than those of pi,2- 

exist singularities higher-order than —1/2. Fig. [T] show 
a CCE whose curvature changes sign, pi, P2 connect 
points with antiparallel velocities and same values of cur- 
vature. Therefore the joint curvature is nonzero accord- 
ing to Eq. (fl~3|) . Therefore they are singularities of order 
— 1/2. P3 is a singularity of order —2/3 since its joint 
curvature is zero. 




FIG. 2. Typical CCE of the d-wave superconductor with 
Bogoliubov quasiparticle energy lower than the spectral gap. 
Four banana-shaped close curves are shown with higher-order 
singularities pi.,3,5,7. 

In d-wave superconductors and iron-based supercon- 
ductors with high anisotropy 0, ||| , some hot spots are 
resolved by FT-STS measurements among all the singu- 
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larities, which are interpreted by the octet model. Some 
of them are obviously close to higher-order singularities. 
Fig. [5] show the four vectors pi, p 3 , p 5 and p 7 of d-wave 
cuprates for illustration, pi 3,5,7 are all higher-order sin- 
gularities because they connect points with parallel ve- 
locities and their joint curvatures are zero. Detailed ex- 
amination shows that the existence of high-order pi,3,5 
are guaranteed by the four-fold symmetry of the CCE 
of d-wave superconductors. On the other hand, the ba- 
nana shape of the CCE, on which the curvature changes 
sign, is necessary for p 7 to be higher ordered, just as P3 
shown in Fig. [TJ Compared with these four, P4 is obvi- 
ously lower-order singular point whose joint curvature is 
nonzero, but the velocities are small so it still has higher 
intensity according to Eq. (TT21 than other singularities 
except Pi, 3,5,7. Similar conclusion can also be applied to 

P2,6- 

For comparison, the octet vector q5 are expected to 
connects the tips of the two banana-shaped CCE in the 
first and fourth quarter. Apparent difference between P5 
and qs can be clearly seen, qs is actually a lower-order 
singularity because the velocities connected by it are an- 
tiparallel resulting in nonzero joint curvature. In fact 
the octet model is an approximate model whose validity 
is guaranteed by the strong isotropy of the energy disper- 
sion, which leads to the fact that the octet vectors are 
close to the higher-order singularities. 
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APPENDIX A 

The CCE governed by Eq. ([3]) are supposed to be com- 
posed of N closed curves, each of which can be param- 
eterized k % x = /i(c) and k % y — 3^(c). i = 1,2, ...,iV the 

curve label and c is the parameter. (c), <j£,(c)) is the 
tangent vector of the i-th curve. Accordingly to Eq. ^ 
we have the NxN one-parameter families of curves, each 
of which is given by 



pi?{t,c)=tiv)-fi(c), 



Py 3 (t,c)=gl(t)~gi(c) 
whose envelope is determined by 



i,j = l,2,...,N (14) 



This indicates that the two tangent vectors are parallel, 
which is equivalent to Eq. (J5J) because the correspond- 
ing group velocities are orthogonal to their corresponding 
tangent vectors. 



APPENDIX B 

To find the divergence behavior approaching to the 
singularities, we transform the above line integral to 
complex contour integral over counterclockwise unit cir- 
cle by letting z = e lt , K u (z) = {F UJ (z) 1 G UJ (z)) = 
(/«(*), 3* (*)) = M*)> V u (z) = v u (t), 



R(p,uj) 



Pu, 



dz 



27ri J\ x \=i zV u (z)[u) - E K ^( z)+p + i0+] ' 



Using the Cauchy residue theorem, we have 



(15) 



zVu,(z)[uj - £ Kw ( z )+ P + i0+] 

1 



\Zk\ < 1}, 

,N < 1}, 



(16) 



where the sum in the first equality is over those k for 
which Zfc is inside the unit circle. In the second equal- 
ity we intentionally separate two parts of contributions. 
One part is denote by C representing the sum of residues 
from poles of [zT4,(z)] _1 which do not lead to singular 
contributions. The other is a summation over the poles 
of [to — H w (p, z) + iQ + ]^ 1 which is infinitely close to the 



zeros of 

w-H w (p,z)=0. (17) 

Eq. (flTf is the complex version of Eq. (fT0| according 
to z = e lt . The real roots of Eq. (fTUJ) correspond to 
the complex zeros of Eq. (IT7|) exactly on the unit circle. 
Complex roots of Eq. (1101) with positive imaginary part 
give rise to complex zeros within the unit circle. Since the 
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complex zeros of Eq. (|17l) correspond to simple poles of 
Eq. (IT6l) . they do not lead to singular contributions. The 
only source of divergence originates from the degenerate 
real roots of Eq. (ITU1) . 



The existence of degenerate real roots demands that 
Eq. (|T0|) and (fTTj) can be simultaneously satisfied, which 
is realized if (p,w) is a singular point. For this case we 
have from Eq. p^|) . 



i2(p,w)~ - 



^Rcs 



2-^k„(t )+p 



V ZZn ' 



\z ± \ < 1 



(18) 



Clear i0 + plays the role of determining only one of the 
two roots contributes to the integral. The result is 



fl(p,w)~- 



«w(*o) 



sgn(£i 



Mt )+pJ 



M*o)+p 



1/2 



(19) 



which indicate that on any singularity the value of R is 
infinite. This result also gives clue on the fact that the 
singularities are the common boundaries for the real and 
imaginary parts to diverge. 



APPENDIX C 



We have applied Eq. (fPOj) and (fTTj) in the derivation. 
From the above equation we find new roots at the vicinity 
of the degenerate root to, 



t± = to ± (-A) 



1/2 



(21) 



where 



Here we give detailed derivation of the divergent be- 
havior of R(p, lo) when any singulary is approached. If 
(p,w) is a singularity, then the three coupled Eqs. (J3|), 
((4]), and ([5]) determine a touching point ko. Or equiva- 
lently Eq. ([TO]) has a degenerate root to satisfying (fTT|) . 
For any give u, when approaching to certain singularity 
p according to R(p + S,u>) where \8\ ~ 0, we need to 
estimate how the degenerate root to varies with 8. To 
leading order of 8 we have the pole equation, 



= uj - K 



k„(t)+p+<5 



u ~ ^k„(t)+p _ <*> ' Vi? k ^(t) +p 



1 ^ 



x ^(t u ) i P (* - W + 8 ■ VK kai(t „) | p 



= 28- V^ (to)+p 



(22) 



Mt )+p 



If A > 0, t± = to±i\&\ 1/2 . The two CCE do not intersect 
at ko = k w (to) any more. This indicates that we are 
approaching the singularity p from the region with less 
number of distinct intersection points. Transforming to 
the complex plane, z + = e lt+ is within the unit circle and 
(20) therefore contribute to the integral according to Eq. (|16l) . 



R(p + 8,co) 



v u (t )zo 

V u (to) 



Res 



Pk„(t o)+P (^f ) 2 + V^ (to)+p ■ 8 - »0+' 
sgn(£ ; klJ (to)+p) 



sgn(-Ek„(t )+ P 



(23) 



|V£ ko ||V£ ko+p | 



2-EK,(t )+ P <5 ' n k + P 



This results indicate that when approaching the singular- real part of R(p + 8, cj), i.c FT-LDOS, diverges according 
ity from the region of p where the real roots of Eqs. © to |<$| -1 / 2 where the imaginary part, i.e. JDOS does not. 
and ([!]) is less than the other side of the singularities the 
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If A < 0, t± = to ± | A I 1 / 2 meaning that the degenerate 
real root to now become two distinct real roots. This 
indicates that we are approaching the singularity p from 
the region with more number of distinct roots. For this 
case the infinitesimal imaginary number i0 + plays the 
role of determining which one of the two poles z± — e lt± 
lies in the unit circle in complex plane. The result is 
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Combining Eqs. (|23| and ([24]) together and making 
use of the result of E given in Appendix D, we obtain 
Eq. CCH). 



APPENDIX D 



Considering that to is a degenerate root satisfying 
Eqs. (Unjl and (fTT]). we have 
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where f> is the joint curvature defined according to 
Eq. (|13l) and n k = r^f|n the unit direction vector of 
the group velocity. Since k and k + p satisfy Eqs. (J3|), 
(|4]) and (JSJ). Therefore n ko+p • n ko = 1 if the group ve- 
locities are parallel and —1 if antiparallel. K k represents 
the curvature at point k on the CCE, 
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(26) 



where E% = fg-, and Ej^ = 

Obviously if q = 
erate root with multiplicity 3, which definitely gives rise 
to higher-order divergence. 



k - -m-dk; with a >P = x >y- 

0, then E = and t will be a degen- 
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